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1.0 INTRODUCTION 


1.1 Background 

The geostationary equatorial orbit (GEO) has been extensively used 
during the past 10-15 years. Satellites in this orbit appear nearly 
stationary above the earth arid are ideal for applications in communica- 
tions and earth observations. 

With the increasing applications of the GEO? the available space in 
the orbital arc is now being depleted. For example? the NASA Space 
Transportation System is expected to launch approximately 160 geosta- 
tionary missions during the next ten years (Reference 1). The satellites 
need to be separated by a few degrees in order to avoid physical and/or 
rad i o frequency interference. 

Because of the extensive planned geostationary mission traffic and 
the limited number of slots available? international organizations are 
attempting to gain authority over slot allocations. The U.S. has no 
inherent claim to geostationary "real estate" since it has no territory- 
over the equator. The geostationary arc of interest to the domestic U.S. 
(55 deg W. longitude to 135 deg W. longitude) may in the future be shared 


with South American countries. Also* certain futuristic projects (such as 


the solar power satellites and space manufacturing) will require large 
amounts of the geosynchronous or geostationary arc. 

For the reasons cited above* it will be important in the future to 
pack more satellites into useful geosynchronous orbits. A promising pro- 
posal that has been put forth is to use geosynchronous orbits whose planes 
are inclined to the equator. A satellite in such an orbit has a ground 
track that is a figure eight. If several satellites are Placed in orbits 
that have the same figure eight ground track* they will all cross the 
equator at the same longitude* called the "gateway". The gateway may 
require only a few degrees of the geostationary arc* while accomodating 
possibly 10-12 satellites in inclined orbits. As the satellites pass 
through the upper (or lower) latitudes* their ground track would appear to 
be nearly circular about a point on the earthy's surface. From this point 
of view, they are called "halo" orbits. 


1.2 Overview of this Report 

Most of the earlier- studies of inclined geosynchronous orbits have 
concerned geometrical relationships with idealised (unperturbed) orbits.. 
In practice* however* the orbits will be perturbed by gravitational and 
non-sravi tat i onal forces. The result is that the orbits will be altered* 
and corrective maneuvers must be executed by on-board control systems. 
Principal orbit perturbations are due to earth oblateness* earth equato- 
rial ellipiticy* sun and moon gravity* and solar radiation pressure. 


For equatorial geosynchronous (geostationary) orbits* the effects of 
these perturbations are well known. However* for inclined geosynchronous 
orbits* less information is known about the orbit perturbations or the 
orbit maintenance delta-V requirements. 

The purpose of this report, therefore, is to provide certain data that 
is needed for preliminary design of inclined geosynchronous missions. In 
part i cu 1 ar * the objectives are to* 

• Study the perturbations on circular orbits with inclinations 
up to 60 degrees. 

• Determine the time histories of the orbital elements for 
periods of up to ten years. 

• Determine the velocity management requirements needed to 
cancel the major perturbing effects. 

• Present a compilation of data on inclined circular 
geosynchronous orbit characteristics. 

Section 1 of this report introduces the problem and gives the scope 
and purose of the study. A summary of the principle results is also 
given. 

Section 2 describes the inertial and earth-fixed (rotating) co- 
ordinate systems* as well as orbit parameters and elements. The complete 
family of geosynchronous orbits is discussed. It is shown that circular, 
inclined geosynchronous orbits comprise only one set in this family. 

Section 3 gives a discussion of the major orbit perturbations, and 
their separate effects on the geosynchronous orbit. 


ORIGINAL, 


Section 4 gives detailed information on the orbit perturbations of 
inclined, circular geosynchronous orbits. The major emphasis is to pro- 
vide time history data of certain orbital elements. 

Section 5 provides a determination of orbit maintenance delta veloc- 
ity requirements to counteract the major orbit perturbations. The purpose 
is to provide order of magnitude estimates, and to show the effects of 
orbit inclination on delta-V. 

Section 6 discusses some of the considerations in mission design for 
a mu 1 1 i sate 1 1 i te system, i.e. a "halo" orbit constellation. 

Section 7 presents some conclusions and recommendations based on 
this study. 

Appendix A gives a bibliography of papers and reports on the subject 
of geosynchronous orbits. There is a wide body of literature on this 
subject. The bibliography includes papers of interest to this study, as 
well as background papers concerned with the fundamentals of orbital 
mechanics . 

Appendix B contains a listing of all the computer programs that were 
developed for this study. These programs are on the Interdata 8/32 com- 
puter at the NASA Johnson Space Center. All routines are throughly 
documented with comment cards. 
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1.3 Applications of Inclined Geosynchronous Orbits 

Several authors have discussed the possible applications of Inclined 
Geosynchronous Orbits (IGO). Generally these orbits are useful for 
applications where there must be a Ions viewing time from the earth. How- 
ever it is not necessary that the satellite be fixed relative to the 
surface of the earth. The following applications have been proposed for 
I GO orbits! 

• Earth observation 

• Solar Power Satellite 

• Communications 

• Military applications 

Applications and around tracks are discussed by Akin (Reference 2>» 
Bielkowicz (Reference 3) and Graf (Reference 4). 


1 . 4 The " Ha 1 o " Orb i t Co n ce pt 

There is an additional application of IGO that has been proposed by 
Fielder (Reference 5). His suggestion is to establish a constellation of 
satellites in inclined circular geosynchronous orbits, such that each 
satellite in the constellation follows the same figure eight ground track 
(see Figure 1). At the higher latitudes (which are of interest for the 
U.S.). the satellites would appear to rotate in a nearly circular orbit 
relative to the earth, hence the term "Halo" orbits. 

The advantage of the "Halo" orbit scheme is that it could signifi- 
cantly increase the number of usable satellites in geosynchronous orbit. 
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Summary of Results 


The general results of this study can be summarised as follows: 

• The general behavior of the geosynchronous orbit per- 
turbations are similar for the inclined (up to 60 deg.) and 
equatorial cases. However? the magnitudes of the perturbing 
accelerations can be different by several orders of magnitude. 

• When the orbit is allowed to have a significant inclination? 

the mission design becomes more complicated. For example? 
assuming the inclination is given? the inertial orientation 
of the orbit —O- must still be specified. Additionally? 
orbit perturbations act in a different way for different 
values of _/"L . 

• The drift in mean longitude is generally less for inclined 
orbits. Likewise? delta-V to cancel this drift will also be 
less for the more inclined orbits. 

• The out-of-plane orbit maintenance delta--V for inclined orbits 
can be as much as four times that for equatorial orbits. Also? 
the amount of delta-V required depends on orbit plane orienta- 
tion relative to inertial space. This means that each satel- 
lite in a "Halo" orbit constellation would have unique 
delta-V requirements. 

• By far the most important source of orbit maintenance delta-V 
is the out-of-plane component. This could be 100 times 
larger than the in-plane delta-V. 
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2.0 ORBIT GEOMETRY AND DEFINITIONS 


In this section* we take the opportunity to review the geometric 
relationships between the geosynchronous orbit and coordinate systems. 
The orbital elements that are used in the following sections will be 
defined here. Also* this section will categoric;-; the various types of 
geos v n c h r o n o u s orbits. 

2.1 Orbital Elements and Coordinate Systems 

Define a rectangular cartesian coordinate system with its origin at 
the center of the earth. The X-axis lies m the equatorial plane and is 
directed toward aries (the vernal equinox). The Z-axis lies along the 
rotational axis of the earth in the direction of the north pole. The 
Y-axis lies in the equatorial plane and completes the right-handed co- 
ordinate system (See Figure 2(A)). 

The following orbital parameters are defined in Figures 2(A) 

Incl i nation* 

Longitude of the ascending node* 

Argument of perigee* 

True anomaly* 

Satellite radial distance from the earth"' s center* 

Bern i -ma.j o r ax i s . 


and 2(B): 

J : 

•XL : 

Ur : 

f : 
' r : 

CL s 
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Apogee distance, 



Perigee distance, 


JE : Eccentric anomaly. 


We also define the orbital eccentricity as 



& CL 


In order to study the satellite motion relative to the rotating 
earth, it is usefu’ to introduce polar coordinates and a rotating co- 
ordinate system that is fixed in the earth. Refer to Figure 2<C). Let 
the X'~axis lie in the equatorial Plane at the intersection of the 
Greenwich meridian and the equator. The Z'-axis and the Z-axis coincide. 
The Y'--axis lies in the equatorial plane and completes the primed 


coordinate system. 

We define the following variables: 

Angle between the X-and X'-axes, generally referred to as 
the "Greenwich Hour Angle", 

The satellite's longitude with respect to the Greenwich 
meridian, 

The satellite's latitude. 

Also, notice that for a satellite in an inclined, circular 24-hour 
orbit, the longitude J\, will have some mean (or average) value over the 
course of a day (Fig. 2(D)). The ground track is a figure eight and will 
cross the equator at a longitude, which is the mean longitude. We refer to 
mean longitude as 7\ ’ a parameter which indicates where the ground track 

is located at anv time. At the instant that the satellite crosses the 


A : 
6 ■ 
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equator soiris north? we have the relation 

L - (p? i~ h. 


1 



Figure 2(a).- Orbit parameters. 
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Figure 2(d).- Satellite groundtrack. 


2.2 Classification of Geosynchronous Orbits 

In this report> the term Geosynchronous Orbit (GO) refers to an orbit 
that has a period of one sideral day? i.e. the satellite's orbital period 
is equal to the period of rotation of the earth relative to inertial 
space. Thus? a satellite in such an orbit appears to follow the same 
ground track? day after day. For these orbits? the inclination arid eccen- 
tricity are not necessarily aero. Therefore? their ground tracks will be 
closed curves of various shapes. Examples of such ground tracks are given 
in References 2? 3? 4? 6 and 7. 

The geosynchronous orbit is the most general type of orbit 
that is under consideration here. Indeed? it defines a family of orbits 
that have the following general characteristics! 

• The ground track is a closed curve. 

• The ground track is repeated every day. 

• The satellite crosses the equator twice each day. 

A subset of the general case includes the circular orbits with an 
arbitrary inclination and longitude of ascending node. These are called 
"Inclined Circular Geosynchronous Orbits" (ICGO>? and are the subject of 
this report. Their ground tracks are the shape of the familiar figure 
eight. 

We suggest the breakdown of the GO family of orbits into sets and 
subsets? as shown in Figure 3. Note that the Equatorial Circular Geo- 
synchronous Orbit (ECGQ) is generally referred to as the Geostationary 
Equatorial Orbit (GEO). A Satellite in such an orbit appears to be 
stationary above the rotating earth. There is only one orbit in this set. 
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Each ECGO satellite is in the same orbits but at a different position 
within the orbit. The satellited position can be specified bv onlv one 
parameter* i.e. its longitude. 

When an orbit is designed to have a. significant inclination and/or 
eccentricity* the number of degrees of freedom for the orbit are 
increased. The degrees of freedom for each GO subset are shown in Table 1 . 
Note that the ECGO set has onlv one degree of freedom* i.e. the mean 
longitude. The ICGO set has two additional degrees of freedom* i.e. the 
the inclination and ascending node. 

The study in this report is concerned with the Inclined Circular 
set of geosynchronous orbits (ICGO). The following analyses in Sections 4 
and 5 assumes that the eccentricity is small and negligible* and that the 
inclination is in the range 0 < I < 60 . 
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FIGURE 3: Breakdown of the GO Family of Orbits 
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3.0 PERTURBING EFFECTS ON GEOSYNCHRONOUS ORBITS 


In this section we give a brief summary of the major orbit perturba- 
tions on geosynchronous satellites. The perterbations of interest are 
those of Ions period that can have an appreciable effect over a period of 
months or years. Short period perturbations (those that have periods 
on the order of a few days) will cause small oscillations about the 
path of the long period trends in the orbital elements. The mathematical 
methods employed specifically exclude these oscillations in order to 
determine the more important Ions period effects. 

The dominant perturbing accelerations on the satellite are: 

• The triaxial earths i.e. earth oblateness and 
equatorial ellipticity. 

• Sun and Moon gravitation. 

• Solar radiation pressure. 

The effects of each perturbation are briefly described in this section. 
Additional background information is given in Reference 6 . The quantita- 
tive effects of these perturbations are described in Section 4. 

3.1 Ellipticity of the Equator 

The earth’' s equatorial ellipticity causes a slow drift (east or west) 

in the satellite's mean longitude . These perturbations are caused 

b v t hi e T 2i t e r m i n t hi e s e o p o t e ri t i a 1 . 

■ — ^ 

Let Y~ be the satellite's position vector in the X?Y ? Z coordinate 
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system in Figure 2(A). Then acceleration due to the seopotential is 

^ IS . 


r - 


cP r 


where 


y - dds. ( 1 - 14). 


The parameter UAe is the earth's gravitational constant. 

Let Kg. be the mean equatorial radius of the earth? and let and 

R. be Legendre functions and associated Legendre functions? respect- 


itvelv. Then the earth's geopotential can be expressed as (Reference 8) 

14-2 T„ (%) n £„ 6 ) + 

+ 4 £ (A- Iip), 

Y\~X ri:zl 

For a triaxial earth? the geopotential becomes 

14 r = 7 X (tyf + 

+ 7 l3 . (Is-'f c**, z(A-'\xi') . 

Observe that the term (oblateness) depends on the satellite's alti- 

tude and latitude. The Tzjl term (equatorial ellipticitv) depends? 
additionally? on the satellite's longitude relative to the earth refer- 
e nee 1 o n s i t u d e * 2X - 

The seopotential parameter 2 i is the angle between the equatorial 
ellipse minor axis and the Greenwich meridian? as shown in Figure 4. The 
non-dimensional coefficients have the values (Reference 8) 

J, = /, OS’^HS 1 X to 1 , X, = ), 77//bX U> . 


I 


Y 



Figure 4.- The Earth's equatorial ellipticity. 
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3.2 Sun* Moon and Earth Oblateness 

The * 2 . term in the geopotential combines with sun and moon gravity 
to cause a long period perturbation of the satellite's orbital plane. The 
period is more than 50 years. 

For I initially zero* the period is 53 years and the amplitude is 
15 deg. This means that after 26 years* the inclination will have 
increased to about 15 deg. 

Section 4.2 describes the motion of inclination for arbitrary initial 
incl i nations. 

3.3 Solar Radiation Pressure 

This non-sravitational perturbation is usually small* but could be 
appreciable for geosynchronous satellites with large solar collectors. 

For the Solar Power Satellite (SPS>* this perturbation has the order of 
magnitude of gravitational terms due to sun* moon and Reference 4). 

The solar radiation pressure acceleration is proport i onal to the satel- 
lite's area-to-we iglvt ratio. The major effects of this perturbation is a 
v e a r 1 y o s c i 1 1 a t i o n of the e c c e n t r i c i t y a n d a r • o t a t ion of the 1 i n e o f 
apsides. Section 4.3 gives orbital element time history data for this 
per tur bat ion . 
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4.0 TIME HISTORIES OF ORBITAL ELEMENTS 

In this section uie study the time histories of the orbital elements 
that are important for geosynchronous missions. We generally relv on 
analytical formulas to generate the time histories. This allows us to 
parameterize a broad spectrum of initial conditions. Also* the formulas 
will be used in Section 5.0 to determine orbit maintenance delta V. 

There are some limitations in the analytical formulas that will be 
used. For example, each set of perturbations discussed in Section 3 is 
taken separately. The coupling (i.e. interaction) of the perturbations are 
ignored since the coupling will be of second order in magnitude compared 
to the principle effects. However, these formulas should provide suffi- 
cient accuracy for mission planning type studies. 

In' Section 4.1 the Ions term effects of earth 1 ' s equatorial ellipti- 
citv ( ’Tj-i- ) are discussed. Section 4.2 discusses the out-of-plane 
perturbations caused by sun, moon and earth oblateness ( U'i. ). Finally, 
Section 4.3 describes the perturbations on the orbit's eccentricity due 
to solar radiation pressure. 
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4.1 Earth-Re 1 ative Mean Longitude 

A oood description of the longitude drift for equatorial saltllites 
is given bv Blitzer in Reference 9. In this section? however? we will use 
the equations from Wagner (Reference 10) since his equations are valid for 
non-zero irtcl i nations. 

Let 7\ be the mean daily longitude of the satellite (Section 3.1). 
This is also the lonsitude of the satellite as it crosses the equator in 
its figure 8 ground track. From Reference 10? we have the acceleration of 
the mean lonsitude siven bv the equation 


9> 


3 


£ Fh, F (^21. ^ ^ 2. 2^) “j 


where 


F; 


2.1 


(s> < 

r z 





1L. ) 


F(i) 


XI 


_ ( 1 + X) 


(4. 1 > 


We define the following physical quantities: 

sLig : Gravitational parameter of the earth? 

: Radius magnitude of the srosvnchronous orbit? 

: Radius of the earth? 

31 j. : Geopotential constant? 

Geo potential constant. 

The numerical values for these quantities are given in Appendix B. 

The function F 0^21. is called the "inclination factor". It has 
the effect of decreasing the acceleration of T\ . Thus? it can be seen 
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from equation (4.1) that the nature of the longitude drift for inclined 
orbits is similar to the equatorial case. 

Notice that only the ellipticitv of the equator ti.e. J^z. term) is 
included in this analysis. All the geopotential terms are given by Wagner 
in equation ( 1 ) of Reference 10. However? we only include the 3*1.2 term 
here because that will give sufficient accuracy for our own purposes? and 
also will demonstrate the nature of the motion. 

In using equation (4.1) to describe the motion of the mean longitude? 
we have made the following assumptions! 

• Terms proportional to the eccentricity can be neglected. That is? 
the orbits are nearly circular. 

• The 7*i term in the geopotential describes the longitude drift 
with sufficient accuracy. 

• There is no coupling of Tn with the sun and moon perturbations. 
This means that the inclination does not change significantly during 
the time of interest. 

The sun and moon have a very small direct effect on the mean longi- 
tude drift. Their direct effect is on the order of the neglected geopo- 
tential terms. However? in the next section it is shown that the inclina- 
tion can change bv several degrees in a few years. Thus? a "mean” (or 
average) inclination should be used in equation (4.1). This should give 
sufficient accuracy for our- purposes here. 

The long-term drift in mean longitude is shown in Figures 5? 6 
and 7. The data presented in these figures was generated from the LDRIFT 
program that is listed in Appendix 8. This program includes the following 
r ou fines! 

LDRIFTs Main program and numerical integration driver. It also prints the 
data to the CRT screen and writes the data to an output file. 

This file is later input to the plot program PLOTIT. 
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L INPUT! Routine for- input and initialisation of constants and common 
blocks. 

LDERIV* Routine to compute the the right side of the differential 
equations <4.1). 

RKF45 : Runge-Kutta numerical integration rout 7 ' . This routine is part 
of the FDS-2 utility library and is not listed in Appendix B» 


DISCUSSION OF FIGURE 5 

(i) These figures show the motion in "phase space"* i.e. the motion of 
longitude drift and drift rate. The arrows indicate the direction 
of the motion. 

<2> Notice that when the drift rate is zero* the drift acceleration is 
not zero. 


(3) Each geosynchronous statellite will be on one of the closed curves 
or "orbits" sho^n in Figure 5. Thus* it is assumed that that a 
satellite will librate and not be in the "circulation" region. 


(4) The optimum location to place and maintain a satellite is at the 

position in phase space where the drift rate is zero. That will be 
on the horizonal line in Figure 5. The satellite will then tend to 
move away in the direction of the arrows. 

<5) There will be similar plots for mean longitudes in the range ISO 
to 360 des. 

<6> The effect of orbital inclination is to squeeze down the the orbits 
in phase ..pace. The total excursion in longitude is not changed* 
but the maximum drift rates will be less for the higher inclina- 
tions. This is represented as decreased orbit correction delta-V 
requirements* as shown in Section 5.1 


DISCUSSION OF FIGURE 6 

(1) These figures show mean longitude as a function of time* for 
various initial longitudes and inclinations. 


(2) Notice that the amplitudes of the oscillations depend very strongly 
on the initial longitude. The stable longitude is near 75 deg. 

Thus the amplitudes in Figures 6(0 and 6(D> are the smallest of 
those shown. For /'p- 1S"° < or loS ° west), the figure would 

show nearly a straight line. 
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(3) The orbit is initialised near the unstable point in Figure 6 < G ) . It 
is seen that the amplitude is very large* but that the motion is 
initially very slow. Thus we can expect verv small fuel require- 
ments in order to maintain the satellite's position near the 
unstable points <-15 and 165 des), 

(4> The effect of an inclination is to reduce the rate of change of the 
mean longitude. 

CUSS I ON OF FIGURE 7 

This figure demonstrates how the period of the libration will 
change as a function of inclination and initial longitude. Notice 
that the function is highly non-linear. Orbits with larger inclin- 
ations have a longer period of libration. 
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FIGURE 6: MEAN LONGITUDE VS. TI 

(A) INITIAL LONGITUDE = 0 DEG. 




















FIGURE 6: MEAN LONGITUDE VS. TIME 

(G) INITIAL LONGITUDE = 165° 





FIGURE 6 ‘ MEAN LONGITUDE VS. TIME 
(h) INITIAL LONGITUDE = 170° 
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INITIAL LONGITUDE RELATIVE TO STABLE POINT (DEG) 
Figure 7 : Libration Period for Inclined Orbits 


-I.,' Motion of the Orb i t Plane 


The differential equations for the motion of the orbit Plane are 
taken from Reference 11 . The equations for inclination X and ascending 


1 1 o de { i are 


(fi /" ^ ; y 


op.SL) 


rA* 


JTL 


•¥ Csv XJTL) # 


( 4 . 2 ) 


i° (t-y CSV 2 l # 

<•„ i E /ct-' • J~ -* 

The s.mbol EL is the eccentric anomaly. For small eccentricities! an 
increase in E of zr corresponds approximately to one dav in time. 

The parameters in the above equations depend on the sravi tational terms of 
the suri! moori! and J2i as well as the obliquity of the ecliptic and the 
altitude of the satellite. The'.- are non-dimensional parameters with the 
val ue r -i 

£ rc , , QQ X i& ) ^ : “ 1 J i~ * # 1 / '■* ^ 7 . 

The derivation of equations (4.2) and the definition of the param- 
eters are given in Reference 11. Also given there is an analytical 
solution that is valid for small inclinations. 

We are interested here in geosynchronous orbits that may have a 
larse inclination. For this case, there is no analytical solution to 
(4.2). However! we can solve them by computing a numerical solution. 

Note that the following integral of motion exists: 


^ ~ o', /Q.A"y . I - p jL 

“ y > 2 Ca>v? LX -/X , 


OK? 
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An advantage of the analysis done in Reference 11 is that concise 
formulas result. The formulas do not provide high precision* but do 
provide enough accuracy to adequately describe the motion. That is our 
purpose in this report* i.e. to determine the nature of the motion of the 
orbital plane. 

In looking at the effects of the perturbations on the orbital plane* 
we first look at the rate of change of X and _/Z . It is of interest to 
compare h i and for different inclinations and nodes. 

It can be shown that for Tp “ O 1 then AT~ •£ > y} deg/' •"ear. 
Compare this with the data in Figure 8. It can be seen that for non-zero 
incl inations* AT is generally less than for the zero inclination case. 
In fact* when XX 0 i&O* * then AT^O . However* this does not 

imply that the orbit will be stable. By observing Figure 9* it is seen 
that A XX can not be zero. Thus as XX moves* it will induce a 
motion in X ■ These considerations will be important in computing 
out-of-plane orbit correction delta-V requirements in Section 5. 

The solutions of (4.2) were obtained on a desk-top calculator* using 
an Euler method for the numerical integration. Time histories for inclina- 
tion are shown in Figures 10 for various initial values of inclination and 
ascending node. 

DISCUSSION OF FIGURE 10 

(1) During the course of 1,0 Years* the inclination can change by up to 
to 8 degrees* either increasing or decreasing. The amount of the 
change depends strongly on the value of the initial ascending node. 
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(2) Only one curve is shown for the X 0 ” O case (Figure .10(A)). This 
is because the ascending node is not defined when the inclination 
is 7.ero. In Reference 11, it is shown that as the inclination 
increases? the node becomes defined at 90 dee east - of the vernal 
esuinox direction and moves west (i.e. — Q_ decreases from 90 deg)* 

(3) The "inclination offset" method for maintaining inclination is 

shown' in Figure 10(A). Notice that when 1° and SLo ~ ^ ? the 

inclination remains less that 1 deg. for a longer period of time 
than when O . 

(4) The general nature of the inclination time histories are similar 
for initial inclinations of 40 and 60 degrees. 

(5) Remember from Section 2 that each satellite orbit in a "figure S 
constellation" has a different ascending node. Thus? we see from 
Figures 10(B) and (C> that the orbital planes of the satellites 

in the constellation will have a complicated motion relative to one 
a n o t h e r . 
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INITIAL ASCENDING NOE 



















4.3 Change in Eccentricity 

As discussed in Section 3.3* the solar radiation pressure will cause 
the orbital eccentricity to change and the line of apsides (once it can be 
defined) to rotate. Ordinarily* the solar radiation pressure has a very 
small effect on the orbit* and is not usually considered in advanced plan- 
ing studies. However* we will study this perturbation here because solar 
radiation pressure could become important for geosynchronous satellites 
with large solar arrays. 

The direct gravitational effects of sun* moon and .J2 on eccentricity 
is small. That is because these terms are proportional to 0 » which is 

considered to be small. Therefore* we include only solar radiation 
pressure effects on orbital eccentr i c i tv. 

As before* we assume that the eccentricity regains rather small 
(less than 0.1). Therefore* we can use the euations in Reference 4. 
Introduce the non-singular elements 

P~ e CtJU? (us-tJTL) J e / °^' 1 • 

These elements remain defined for circular orbits. 

In developing the acceleration of the satellite due to solar radia- 
tion pressure* the following simplifications are made: 

• The satellite can be modeled as a flat surface which is 
perpendicular to the satellite-sun line. 

• The satellite is a 10% reflecting body. 

• The earth moves on a circular orbit about the sun. 

• The solar energy flux is constant in the vicinity of the 
earth . 
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Further details on the above described model are siven in References 4» 
12 and 13. 


After makine the assumptions and simplification described above* 


the differential equations for p and are found in Reference 4 to be 

Jp 


■IE ^ 


£ 


[n.‘ £ q + g Q (% £ -t- 



'& — 


-h 


-[riC ri $ £ <£,?£ (?£+$£)■ 

P 9 (gltQp) (l-Ctf'ff 1 )) (2- (£‘+t? 1 )') /l J * 


(4.5 


f £ ~ I f>$ ( 2 v + 

- jV Cf £ Q - ri‘ 1 0 Q (2 - (£ Vc? 0) Vt f C p £ (t£+Q p) 

- ?p(z£ iQ ^ (1- (q'+£ ‘)) (2 ~ ('£ Xf Q 1 )) A J ^ v ) ‘ 

The above equations make use of the follouins abbreviations: 

£ ~ (l ~ Co? I c^p „_0_ j Q~ — )/l - E^?jp /2-Ao J1 j 

2 -(p l t£ ; ) 


= £ 

SbS.lZ' 


^ Qq . 


Also* define the f 0 1 1 o wine parameters: 



anele of the sun in the ecliptic plane* measured from the 
vernal equinox. 



/i 


(o~h. sr°) 
( S.oL ', X /J y ) 


(4.6) 


is the cross-sectional area of the satellite in square meters? and 
\^/ i s the weight of the satellite in kilograms. The parameters ^5 
and have been previously defined. 

In order to determine the Ions term change in eccentricity? the diff- 
erential equations (4.5) were solved by using a Runge-Kutta numerical int- 
egration method. Element time histories are shown in Figures il?12 and 13. 
The data presented in these figures was generated from the EDRIFT program 
that is listed in Appendix B. This program includes the following 
routines! 

EDRIFT: Main program and numerical integration driver. It also prints the 
data to the CRT screen and writes the data to an output file. 

This file is later input to the plot program PLOTIT. 

EINF'IJT: Routine for input and initialization of constants and common 
bloc k s . 

EDERIV: Routine to compute the right sides of the differential 
equations (4.5). 

RKF45: Runge-Kutta numerical integration routine. 

The data in Figures 11? 12 and 13 was generated based on an area to weight 

. . \-l 


ratio of 1.73 ft) 


for . 


This is an assumed value for the Solar Power 


Satellite (Reference 4). However? since the perturbing parameter 


eter 


depends directly on 


(see equation (4.6))? we can expect that 
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satellites with different h/W will have proportionate behavior as 


shown in Figures lli 12 and 13. 

DISCUSSION OF FIGURE 11 

(1) The purpose of these eight figures is to show the eccentricity time 
history* and how it depends on X and _fL 

(2) When the orbits lie near the ecliptic plane ( XL ^ O >, the inclina- 
ion does not have much effect on eccentricity time history. 

(3) All curves have a period of one year. However, the maximum value of Q. 
can depend strongly on inclination and node. 


DISCUSSION OF FIGURE 12 

(1) This figure can be thought of as a plot of f* and £ in polar 

coordinates (see equations (4.4)). As time increases, the closed 
curve will be traced in the direction of the arrows, starting from P~Oj 

The value of the eccentricity will be the distance from the 
origin. The angle Ur + -O. will be the angle between the radius vector- 
arid the p — o axis. For example, in Figure 12(A). us+-TL begins at 
aero degrees and asymptotically approaches ISO deg., to start the 
cycle again. 

(2) For different inclinations and nodes, the shape of the closed curve 
will change. In the 1=0 case, the curve is slightly flattened. This 
is because the orbit is inclined to the ecliptic plane (the Plane of 
the sun). 

(3) For a node of 90 deg. and inclinations of 20. 40 and 60 deg., the 
curve is flattened. Compare with Figure 11(0. 


DISCUSSION OF FIGURE 13 

The purpose of the figure is to show the 'long term (10 years) motion 
of eccentricity. Note that the motion is periodic with a period of 
one year. Thus, it is only necessary to look at the motion within 
one year. 


52 



Time History 
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Figure 13: Long Term Motion of Eccentricity 



5.0 VELOCITY REQUIREMENTS FOR ORBIT MAINTENANCE 


In this section u»e determine the delta velocity increments that are 
required to offset the orbit perturbations. 

5.1 In-Plane Delta-V 


5.1.1 Mathematical Developments 

This analysis presented in this section makes use of equations and 
results from References 10? and 14. The following two assumptions are 
made! 

(1) The orbit is nearly circular. Thus the delta-V ( AV ) in the 
radial direction is negligible as compared to the tangential 
component. 


(2) There is no coupling between the in-plane and out-of-plane pertur- 
bations. This is a limitation of the analytical perturbation 
method? and was discussed in the previous section. 

The orbit-averaged longitude drift of an inclined circular orbit 
was given in Section 4.1 as . 


(5. 1 > 


where 

Since we 
i n g n o d e 
equat i. on 




~ ^ f-S/W. day) 

are interested in the deviation of the mean longitude of ascend- 

of geosynchronous orbits? it is convenient to linearize the above 

in terms of the variation h ~h defined by the relation 
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/ p. o ) 

\ Sm* • X. / 


/\ “ t 2) /) 

where is the desired mean 1 ona 1 tude of ascending node of the orbit. 


By substitution of ( 5 . 2 ) into ( 5 . 1 ) we get 

tO (* 


/l ~~X "f £“#2 A Z 2. C*&oSl(~?\ 5 "All. *2. (2$ ^ 22^, (5. 3) 

The solution to the corresponding homogeneous equation of ( 5 . 3 ) is 

A U ^ ~h & 2, /<2 X'*? LAzJA j 


where 


U/~ 2r £-2 2. ( x s ~ A 


\~1 * J 2- 


2-i 


The particular integral can be found bv assuming 


A Ap - ^3 

S u b s t i t u t i n g i n 1 0 e a u a t i 0 n ( 5 . 3 ) a n d s 0 1 v i n g f 0 1 - C x we g e t 

C 3 = - T 

T h e r e f o r e 7 t h e c 0 m p 1 e t e s 0 1 u t i 0 ri i s 

f\ 7\ ~ Cf -t C \ /<2A^x CrUjt — “tT 2* (A*r • 

The i n teg r a t i 0 n c 0 n s t a n t s c a n b e f 0 1 j n d b v a s s 1 j rn 1 n s t h e r e i s n 0 
initial orbit injection error? so that i ^ ^ 
and A ?l ~ & at J?t ~ O ' 

( 1 ) For jpj'^o and i\~X ~0 ? then 

C ! - A Jtas” 2 £/\ 3 “ /Lz ) ' 


( 2 ) For — 


and 


^ ^ O » then 


C 2 . — o 


T h e r e f o r e t h e s 0 1 u t i 0 n b e c 0 rn e s ? 


A A ~ (7\s ^\ctx&u/Jr- ^L~J 


( 5 . 4 ) 


Since lA/“ is sr&& 1 1 * equation (5.4) can be further simplified by expanding 
Creew J and inserting t h e expression for lAf . The result is 


47 1 = -/Ij.iXf 1 (%- 7*4 


<5.5> 


In Reference 10 it is shown that the rate of chan ye of the circular 


synchronous orbit radius satisfies the equation 

, /24^ } (~X~ ) 


(5.6) 


where 


r - 

m (jf-) 7,. F(l)n. 

Equation (5.6) can linearised in terms of variations /j and /] ]T , which 
is defined bv the relation 

Y- n + 4 r. 

8 u b s t i t ’ t e t h e e x p r • e s s i o n s i n v o 1 v i. n y Ar arid n into eauat i on ( 5. 6 ) to yet 

A X ~ B xx /C ^ % ( 7 \ 2B Ml c^z( 7 i 6 -^^)A 2 . 0.7) 

By substituting equation (5.51 into (5.7) 

A X — & 2 2 ^ X %2^) ^ 

■j-2 6*2. (7\ b ~ 7 w) [~A i vf -2. (X~ 7 n )y« 


(5.8) 


Solving the differential equation (5.8) 

^ IT 3*2. /OASI i J2. ( /I ^ "/i T. 




A s b e f o r e ? assu m e t h a t t h e r e i s n o i n i t i a 1 i n -j e c t i o n e r r o r ? i . e . f o r 
and A Y'--' 0 ’ then C. ~ O 

We thus have the expression! \T 2 2 , \~1 

sf 2.2.J, 1~^A2 iA C#&(T\< > ''~A 3.Z.)/ (5.10) 


4 ^ — ^72 
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V 


To determine the station-keeping requirements? consider a 

satellite which has been drifting? such that A~A = &~X 0 ■ by 

using equation <5. 5)? the time of drift can be found. At that time it 
will have a radial distance of + 4r ■ For a near circular orbit the 


velocity relative to the inertial space is given bv 

Vf - (ts y) 4><s ~ Xs ) <5.ii> 

where A "F is the variation in true anomaly along the orbit? and 
is the rotational rate of the earth. 

To cause the satellite to return to the original value of 5 ? 

it will be necessary to transfer to a circular orbit of radius ~~ A K 

with a velocity relative to the inertial space given by 

^ - (Y'.-A'r) <4 + r s A f . < 5 . 12 , 

In this new orbit? the satellite will drift back to ? and then begin 


(5. 12) 


the cycle again. (Note that this takes advantage of the linearity for 
for small deviations.) 

The minimum energy coplanar circular orbits transfer is a Hohmann 
transfer ellipse with apogee and perigee distances of % s i-AY' and » 

respectively. The two velocity impulses are 




1+ MAL 

Ks-AY- 


A 1/ = 


V, AL . 
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Vl:: V 

Jit. 


v 1 -b/ ri+or 

' Ks-AY 


-\V 2 ' 


(5. 14) 


Note that these two impulses are applied at one half a sidereal day apart. 
By substituting equations (5.11) and (5.12) into (5.1'3)^arid (5.14>» we set 

A !/, = [A)r<p c + -jrfaTfc -AfAr] } 

av % -\ f 4h ic ~ \ A f ^ AY ^‘ 

Lmn ~ 

By neglecting second and higher order terms of the small variations* we 

9et A v, = - -r a r <P? j 

#. 

A K = - i A K ^ • 

The station-keeping hV requirement is then t 

AV = ) A bj ) i- j AYij = J At fi e j . <s.i5) 

The time period for the cycle is equal to T * where from (5.5) 

T=\a A»/(- 

In equation (5.10) the second term in the right hand side is on 

, -J 

the order of jO of the first term. If we the neglect second term in 
equation (5.10>* we can write V as . 

' Mr = l± 0 ai ^2 <Pe I j <5.i4> 


X /<2 ^ 


which is independent of the time period of the correction cycle. 
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Numerical Results 


5. 1.2 

The equations to compute in-plane delta V have been implemented in 
program DELTAV (Appendix B). Results from that program are given here. 

The irr-Plane delta V requirements are shown in Figure 14. The longi- 
tude is measured relative to the stable points <105 deg W and 75 deg E). 
There will be four longitudes which require zero delta V (the two stable 
and two unstable equilibrium positions). The overall effect for inclined 

orbits is to decrease the delta V requirements. This goes back to the 

• # 

"inclination factor" in the expression for . 








5 . 2 0 u t — c< i — P 1 a. n e De 1 t a V 


This section will discuss the delta V requirements to correct for 
the perturbations of the orbital plane. These perturbations were discuss- 
ed in Section 4.2. 

5.2.1 Mathematical Developments 

The normal to the orbit plane can be represented in inertial coordi- 


nates as 


•A A r * 

— As JT. — £ /oX^ X XL ~hXk CooT 


* A p 

where As i and Afc are unit vectors in the XiYiZ axes? respectively. 

If the orbital plane has been perturbed through the small ans 1 e sJtXL and 

AT i then the normal to this perturbed plane can be represented in 
i n e r t i a 1 c o o r d l n a t e s a s 

U-2^( T A^'^(-' Sl+AS ^ Xo^(sttJ-n.) + 

+Jk Cxt at) . 

Expa.ndinsi a. rid neg lectins higher order terms in AX a n d A XL : 

u'~ sC (fri™ T ah™ xl -t Axl&X^T c^xl -t-AI c+aT/oXoXi) t* 

-b ^ (“XUssjX C^aXl. ~t A XL /OX^T a2ass> XL — AT C**XX XL) ~b 

~b Jft [c^&X — AX ssXnT') , 

T h e a n g u 1 a r rn o rn e n t u rri o f a s e o s v n c h r o n o u s s a t e I 1 i t e i n t h e o r i g i n a 1 
orbital plane is — > * A 

Jl = r f y . 

The angular momentum of a geosynchronous satellite in the perturbed plane 

is = b-P of'. 

Therefore the change in the angular momentum vector is 

^ jr -p j As (A XI /2X^X C*s 2 XL -j- AT C*aT /oX^ jQ.') ~h 
A XL /XssTX &Xs) s/I — AT C*eT C^aXL) 


-A 

~ Xti. 
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The magnitude of ^ 5 can be approxirnted bv — i Vl 

jA^I = h l f [(AT- A-fL a^lY+2 AT ASia^IJ . 

The relation betiueen A-Al and A V is given bv 

A Jb = r-A V. 

If we make the following assumptions « 

AI=IAT , AJi = ~n- AT ) 

then the delta V requirement per unit time is given bv 

j~ - r f i') -h 2 J-/i 

.# * 

where J[ and _jQ. are given bv equation (4.2). 
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Numerical Results 


The equations to compute out-of-pl ane delta V requirements have been 
implemented in program DELTAVN (Appendix B). Results from that program 
are given here. 

The out-of --plane delta. V requirements are presented graphically in 
Figures li'(A) and lST(B). In the first case, the delta V is given as a 
function of inclination. In the second case, delta V is given as a 
f u n c t i o n o f a s c e n d i n s n o d e „ 

DISCUSSION OF FIGURE 1 <T 

( 1 ) It in seen that the delta, V requ i rernerit s for our -of -pi ane orbit 
corrections are nearly 100 times the maximum delta V of in-Plane 
corrections. Therefore, these results clearly predominate when consi' 
dering mission feasibility. 

( 2 ) When 1 = 0 , there is no distinction between A V for different values 
of JX • This is because _ 2 TL becomes undefined. (Figure li’(A)) 

( 3 ) In Figure iSTA), the curve for JX-o s h ci u 1 d i n t e r s e c t t h e a b c i s s a 
at 1 = 7.3 deg. It fails to do that because of the granularity of the 
data. Remember that for 1 = 7.3 deg. and — CU^O , the orbit Plane is 
in equilibrium, and no delta V is required. 


(4) Inclined orbits require much more delta V than do the nominal lv 
equitorial case. 

(5) For a nominal inclination of 40 des. (desirable for "Halo" orbits)? 

the delta V requirements are near 1 y the same for all _/L- . However 

these orbits require nearly the maximum delta V. 


75 



FIGURE 15: OUT-OF-PLANE VELOCITY CORRECTION 






FIGURE 15: OUT-OF-PLANE VELOCITY CORRECTION 

(B) DELTA V VS ASCENDING NODE 





6.0 DISCUSSION OF MULT I SATELLITE SYSTEMS 


6.1 Problem Areas 

It has been proposed that Inclined Circular Geosynchronous 
<Fisure 8) orbits may be useful as a means to significantly increase 
the number of 24-hr. satellites in space. Several authors have investi 
gated the satellite-ground station geometry for such orbits. 

This report has been concerned with the orbital mechanics of 
inclined? circular geosynchronous orbits. It was found that station- 
keeping velocity requirements could be an order of magnitude larger 
•than for equatorial orbits. 

Whereas previous studies on this topic have been concerned with 
the utility of a single satellite in a nonstat i onarv geosynchronous 
orbit? additional studies need to be concerned with a system of such 
satellites. The emphasis should be on mission design concepts that 
illustrate the practicality of this approach. 

We suggest an extension of the work on application* of nonsta— 
tationary geosynchronous (24— hr. > orbits? with the following primary 
objectives: 

(1) Study the applications of geosynchronous orbits that have 
arbitrary eccentricity and inclination. 

(2) Develop mission design concepts for multi-satellite 
"constellations" in geosynchronous orbits. 
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Mission Design Concepts 


1 


An appropriate next step in this study would be to develop mission 
design concepts for a "constellation" of satellites. The following are 
some suggested study areas! 

<i> Develop mission design objectives and groundrules for a 
constellation of satellites in geosynchronous orbits. 

Take into consideration! 

• Individual user requirements. 

• Ground station tracking requirements. 

• Launch vehicle capabilities. 

( 2> Determine the number of satellites required in order that a 

given ground location has continuous coverage. Consideration 
must be given to optimum eccentricity and inclination of the 
orbit. This part of the mission design should consider the 
longitude and latitude of the ground stations, as well as the 
equatorial crossing longitude of the satellites. 

(3) Develop an "inclination-offset" technique that minimises out-of- 
plane orbit correction requirements. This concept has been used 
with existing equatorial satellites, and must be generalised to 
the case of inclined orbits. There is an additional degree of /. 

freedom to this problem since each satellite in a system has a 
different equatorial crossing point relative to inertial space. i 

<4> Study the relative motion between satellites and determine the 

delta V requirements to maintain desired orbit positions. ! 

(5) Investigate the problem of collision and interference avoidance i 

at the equatorial crossing. 
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